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Abstract 

We apply the nested algebraic Bethe ansatz method to solve the eigenvalue problem for the SU(4) 
extension of the Hubbard model. The Hamiltonian is equivalent to the SU(4) graded permutation 
operator. The graded Yang-Baxter equation and the graded Quantum Inverse Scattering Method are 
used to obtain the eigenvalue of the SU(4) extension of the Hubbard model. 
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1 Introduction 

Since the discovery of high-temperature superconductivity, much more attention has been paied to the 
theoretical mechanism for such phenomena. Most proposals concern about the Hubbard model and the 
t — J model [|l|. C.N.Yang |^ advocated the importance of ?7-pairing mechanism and the property of 
off-diagonal long-range order (ODLRO) for the eigenfunctions in superconductivity. And it is stated that 
the wave function of BCS theory do have the ODLRO. Essler, Korepin and Schoutens proposed an 
SU {2\2) extended Hubbard model with the 77-pairing symmetry. They showed that the eigenstates of this 
extended Hubbard model exhibit ODLRO and is thus superconducting. A more general extension of the 
Hubbard model with 77-pairing is also proposed in ref.Q, see ref.[^ for some related works. One of the 
authors showed that the SU{2\2) Hamiltonian can be constructed by the coupled SU{1\1) chains, and 
proposed a new S'J7(4) extension of the Hubbard model It can be proved that the eigenstatc which 
shows ODLRO of the SU{2\2) case is also an eigenstatc for the SU{4:) extension of the Hubbard model. 
Thus we present another extension of the Hubbard model, its eigenstatc possesses the ODLRO property. 
The SU{4:) extension of the Hubbard model is constructed by two coupled SU{2) chains. To some extent, 
it is connected with the spin-ladders which are introduced to describe the quasi-one-dimensional materials 
in the superconductivity 0. Recently, a series of exact solvable spin-ladders including coupled t—J model 
were proposed^ . Different from the SU (4) spin-ladder model, the SU (4) extension of the Hubbard model 
has the Hubbard-like interactions instead of the rung interactions in the ladder models. And another 
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motivation for us to study the SU (4) extension of the Hubbard model is that it is the fermion version of 
the SU{4) spin chain with orbital degeneracy|^, p^ . 

In this paper, we shall give a detailed calculation for the eigenvalues of the S'C/(4) extension of the 
Hubbard model. This extension of the Hubbard model is a fermion chain, so we naturally introduce the 
graded Quantum Inverse Scattering Method (QISM) to solve the eigenvalue problem. The Hamiltonian 
is presented to be an SU{4) graded permutation operator. We thus start from the rational SU{A) R- 
matrixt. By using the graded method, we show that the Hamiltonian can be obtained from the transfer 
matrix with periodic boundary conditions. Then, we can apply the graded nested Bethe ansatz method 
to find exactly the eigenvalues of the Hamiltonain. 

The paper is organized as follows: We shall introduce our model and fix notations in the next section. 
We then shall prove that the one-dimensional Hamiltonian is equal to the graded permutation operator 
and the system is integrable in section 3. In section 4, the graded algebraic Bethe ansatz method will 
be used to obtain the eigenvalue of the Hamiltonian. the corresponding Bethe ansatz equations will be 
presented. Finally, we give a brief analysis of the Bethe ansatz equations in the last section. 



2 Description of the Model 

The Hamiltonian of the SU{A) extension of the Hubbard model is written as: 

^ 1 1 

H = -H^ + uY,{n,, ~ -){n,,, - -I (1) 

where — J2jk ^jk coupled SU (2) fermion chains, and the summation takes for nearest- 

neighbour sites. The local Hamiltonian of the coupled SU{2) fermion chain has the form 

X [c]^Cki + cl^Cji + 2njinki - n^i - rifej, + 1]. (2) 

The operator cja and cj^ represent the annihilation and creation operators of electrons with spin cr =1, ], 
on a lattice site j. These operators are canonical Fermion operators satisfying anticommutation relations 
expressed as 

{cL' Cjr} = S.jSar, {c^„ : Cjr} = {cj^, c]^} = 0. (3) 

We denote by rij^ = c}-^Cja the number operator for the electrons on a site j with spin a. 
Explicitly, the local Hamiltonian (j|) can be written as 

^% = 51 [(c]<TCfc^+4<TCj'T)(l-'^j,-a)(l-nfe,_<,)-t-(c]^Cfe<,-f 4^Cj^)nj-_^nfc,_^] 
-f - njx)(nfc| - nki) + ^(nj - l){nk - 1) 
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where we denote by nj = X)<T=t I ^^"^ number operator for the electrons on a site j. The Hamihonian 
is invariant under spin-reflection Cj| <-> Cji- But it does not have the invariant property under 
particle-hole replacement cj^ <-> Cj^. There are four kinds of state at a given site: 

two of them are fermionic and the other two are bosonic. The state | ti> represents that an electron-pair 
is localized on a single lattice site. The formation of such pairs, called localons, is considered to be a 
mechanism to create 'Cooper pairs']^. 

We introduce here some generators. The spin operators S = X^jli ^j^ ^'^ ^'^^ ^ 

Sj=c]^cji, S]=c]^Cj^, S^j =\{nj^ -riji), (6) 

form an SU{2) algebra [S*, S'^] — 2S^, S^] — , [S, S^] = —S. Since the fcrmion operator is grassmann 
odd, we find that the above spin operators are grassmann even (bosonic). The yy-pairing generators also 
form an SU{2) algebra [r|,r|^ = 2ri^, [rj^rj^] — 77^, [r;,??^] — —77 with 

Vj=Cj^cji, 77] =c]^c]^, r;| = -inj + i. (7) 

We see that those operators be grassmann even. The generator Xj — (rij-f — i ) {uji — ^ ) is also a grassmann 
even operator. We further introduce eight grassmann odd generators. They take the following form: 

Qj<y ~ ^'.-o-Cjcr, = "-j,-o-c]o-, (8) 

with a =1, [ representing spin up and spin down respectively. We shall focus in this paper on the one- 
dimensional case with periodic boundary conditions, that means the nearest-neighbour sites of site j are 
J ± 1. In terms of these generators, the one-dimensional Hamiltonian ff^j+i is written as 

+ + s^l'yl+i - s]s,+, - s,s]+, + 2s^sj+, + ax,x,+, + i. (9) 

Next, we shall show that the Hamiltonian -ffjj+i is equal to the graded SU{A) permutation operator. 
Because it is the graded version of the permutation operator, note that this Hamiltonian does not preserve 
the property of SU{A) invariant, for example, we do not have relation [-ff*^, S] = 0. 

We can use 4x4 matrices to express the above generators. The represention can be realized as 

bj - , - , - -{Ej -Ej ), 

Qn = -E^\ Q}^^-Ej\ Q,i = Ef, Q]^ = Ef, (10) 
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where the matix E" acts on the j-th space with its elements defined as {E"^ = Sa^Sp^. We remark 
that the matrix E"^^ is not a conventional matrix but a supermatrix with grassmann number + ^f}^ 
where = represents grassmann even (boson), and = 1 represents grassmann odd (fermion). From 
the representations (|l0[), we find that we should choose the grading ei = 62 = 1, 63 = 64 = so that the 
representations have correct grassmann number. 

Considering the above representation, we find that the Hamiltonian can be written as 

- ^(-1)'"^^ ®^^i+E^r®^^Ti- (11) 

The right hand side the this relation is just the graded SU (4) permutation operator. Now, we define the 
super (graded) tensor-product as 

[A^BUjm = {-1)^'-^+'"^'^ A^nkBji. (12) 

According to this definition, the non-zero elements of the local Hamiltonian read as: 

rrrO 17"^ _ / (-l)'°'''i5ai^'5/3^, « 7^ /?, /-, -^A 

y^n+iUp - a = l3 = -i = v. ^^'^> 

Thus we have proved that the Hamiltonian is equal to the graded permutation operator, and 
because the Hubbard interaction term C^X^j -^j commute with , they have common eigenvectors. We 
then can use graded algebraic Bethe ansatz method to find the eigenvalue of the Hamiltonian (|^). We 
point out that the Hamiltonian is neither the usual S'C/(4) chain nor the super SU {2\2) chain. 

3 The Integrability of the Model 

We begin with the SU (4) rational R-matrix. The non-zero elements of the matrix have the form 

u + i, 

h a^h, (14) 

where indices a, h take values 1, • • • , 4, i = T is the cross parameter. We know this R-matrix satisfies 
the usual Yang-Baxter equation 

-Rl2(wi - "2)^13(^1 - W3)^23("2 - U3) = ^23("2 " U3.)Ri'i[ui - 1*3)^12(^1 - U2) . (15) 

As in section 2, we still choose the grading ei = £2 = 1, — — 0, that means the grading is FFBB. 
Introducing a diagonal matrix /^J? = {—lY'^'^'^ab^cd, we modify the original R-matrix to the following 
form 

R{u) = IR{u). (16) 

For the non-zero elements of the R-matrix i?^^, we have a relation + + Ec + = 0. Thus we find 
the new defined R-matrix satisfies the graded Yang-Baxter equation which reads explicitely as 

R{n-vt\'^^R{u)l^^:R{v)llll{^^^^^^ (17) 
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We remark that, we have already graded the R-matrix, that means the matrix is now a supermatrix with 
grassmann numbers. For u = 0, i?i2(0) = iPi2, and the elements of P12 read [Pi2]at = {^^y°''^''Sad^bc- 
We know that this P12 is the super permuation operator corresponding to group SU{2\2) with FFBB 
grading. 

In the framework of the QISM, we can construct the L operator from the R-matrix as: 

Laq{u) = Raqiu), (18) 

where a represents the auxiliary space and q represents the quantum space. Thus we have the graded 
Yang-Baxter relation 

i?i2(w - v)Li{u)L2iv) = L2{v)Li{u)Ri2{u - v). (19) 

Here the tensor product is understood to be in the sense of super tensor product defined in relation (p^. 
In the rest of this paper, all tensor products are in the super sense. 

Following the standard method of QISM, the row-to- row monodromy matrix Tn{u) is defined as the 
matrix product over the N operators on all sites of the lattice, 

Ta{u) = LaN{u)LaN~l{u) ' ' ' Lal{u), (20) 

where a still represents the auxiliary space. Explicitely we write 

By repeatedly using the Yang-Baxter relation (p^), one can prove that the monodromy matrix also 
satisfies the Yang-Baxter relation 

Ri2{u - v)T,{u)T2{v) = T2{v)T,{u)R,2{u - v). (22) 

We consider in this paper the periodic boundary condition. The transfer matrix t{u) of this model is 
defined as the supertrace of the monodromy matrix in the auxiliary space. In general case, the supertrace 
is defined as 

tiu) = strT{u) = Y,{-^Y^T{u)aa- (23) 

As a consequence of the Yang-Baxter relation ( |2^ ) and the unitarity property of the R-matrix 

Ri2{u)R2i{-u) = (i + u){i - u) ■ id., (24) 

we can prove that the transfer matrix commutes each other for different spectral parameters. 

[i(u),t(«)] =0 (25) 

Generally in this sense the model is integrable. Expanding the transfer matrix in the powers of u, we can 
find conserved quantites, the first non-trivial conserved equantity is the Hamiltonian. 
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For the R-matrix under consideration, the Hamiltonian can be obtained by taking the logarithmic 
derivative of the transfer matrix at the zero spectral parameter, 



ii',+i(0). 



(26) 



Explicitly, we write the local Hamiltonian H',:,-^ obtained above in a matrix form 
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77124 



^31 



p41 
77142 

^1 
77144 



\ 



(27) 



where iJJ^i acts on j + 1-th space. We can find this Hamiltonian is exactly the Hamiltonian H^j+i (p^. 
We show explicitely here what we mentioned above, this Hamiltonian is not the usual SU{4) Hamiltonian 
because we have negative signs in the above relation, it is also different from the Hamiltonian of the 
supergroup SU{2\2) for all diagonal elements are positive here. We have proved the Hamiltonian under 
consideration can be obtained from the transfer matrix of the exactly solvable model and therefore is 
integrable. 

4 Graded Algebraic Bethe Ansatz Method 



We use the nested algebraic Bcthc ansatz method to find the eigenvalues of the transfer matrix. We start 
from the graded rational SU (4) R-matrix with FFBB grading. Explicit form of the R-matrix is 



R{u) = 
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(28) 



This R-matrix satisfies the graded Yang-Baxter equation (p^. In the framework of the QISM, the 
L-operator is defined by t 



( u- h{u)E]^ 



iEf 



-lEf 
- b{u)E 

lEf 

ieI^ 



22 

3 



lEf 
lEf 
u + lE^'^ 
lEf 



iEf 
iEf 
iEf 
^^Ef 



\ 



(29) 
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where E"^ acts on the j-th quantum space, and we use notation b{u) = 2u + i. Considering the repre- 
sentations of the operator E"^ (|l^), we can also write L-operator in the following form 



Lj{u) = 
/ u-b{u)il 



-iS] 



We can also write it as 



4 ^3 



iQjl 

SI) iQ^^ 

u + i{\ + Xj+f]l) 



ir]} 



(30) 



Lj{u) = 

I u-b{u){nji~nj^nji) 



Hl-"iT)S-x 



«(l-"ii)c]| 



i{l-nj^)cji 

*(l-'^ji)cjT 
u + i{l — Tij + nj^riji) 



(31) 



We choose the local vacuum state as \vac >j= (0, 0, 0, 1)*. Actually this local vacuum state is equal to 
the state | Ti>i- Acting the L-operator on this local vacuum state, we have 



Lj{u)\vac >j = 



( 

V 

/ "I n>, 






u 



-^E]^ 



^Ef 



\ 

u 

M 

?f u + i ) 



u\ u>, 



\vac >j 



(32) 



z|0>, {u + i)\n>i J 



Define the vacuum state as Ivac >- 



This vacuum state is the full-filled state just like the 



'Dirac sea'. Using the standard QISM, we denote the monodromy matrix defined in (EG) as follows. 



Tin) 



( (u) Ai2{u) Ai3(u) Bi{u) \ 

A2l{u) A22{u) A25{u) B2{u) 

A^liu) A32(u) A33(u) S3(m) 

V Ci(u) C2{u) C3(u) D{u) j 



(33) 



The transfer matrix with periodic boundary conditions is thus written explicitely as 

t(u) = D{U) + A33(u) - A22{U) - An(u) = D{U) + V(-l)'''^aa(u). 



(34) 



We write the last equation for convenience in using the nested Bethe ansatz method later. With the help 
of the definition of the monodromy matrix (pO) and the result in relation (B2h , we find that the action of 
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the monodromy matrix on the vacuum state is 
T(u)\vac >= 



/ 







































{u + 



\N 



\vac > 



(35) 



From the Yang-Baxter relation (p2|), we can find the commutation relations which are necessary for the 
algebraic Bethe ansatz method, 



Diu)C,iv) 



?! — ?/ -I- 7 1 

—Cc{v)D{u) Cc{u)D{v), 



V — u 
be 



Aab{u)CM - V ,b,i 



U — V 



u — V 



(36) 



-Cb{u)Aac{v), (37) 



rW(v - 



(38) 



V — U + I 

Here all indices take values 1,2,3, and matix i?'^^^ is equal to the original R-matrix ( p8| ) with indices being 
just 1,2,3. We assume the eigenvectors of the transfer matrix as 



(39) 



where F'^^'"'^"^ is a function of the spectral parameters Vj. Operating the transfer matrix ( p^ ) on this 
assumed eigenvector, applying repeatedly the commutation relations and using the result (35), we have 

t{u)C,M)Cc.Av2) ■ ■■C,Jv^)\vac > F^--'", 

= [D{u) + i~lY-Aaa{u)]CcAvi)CcAv2) ' ■ ■ Cc^{Vra)\vaC > F^'-'-^ 



Vj ~ U ~\~ I 



C,Avi)C,Av2) ■ --Cc^Mlvac > F'^- 



+u.t., 



(40) 



where u.t. means unwanted terms, and t*^^) is the nested transfer matrix. The eigenvalue problem of t(u) 
becomes now the eigenvalue problem of t'-^K Denote the eigenvalues of t and A^'' by A and A*^^^. In order 
to cancel the unwanted terms, we need the following Bethe ansatz equations 

m 

Y[{vk-v,-i){vk+i)^ +vj^A^^\vk)^0, fc = l,2,---,m. (41) 
i=i 

The nested transfer matrix is written explicitely as 

^_2yea+ebi)edi+(«Q+e62)<^<i2'l ^-(«a+e6„_l )ed„,_i _ ^42) 

For the non-zero elements of the R^^\uY^, we still have f.a + (-b + (-c + ^d — 0, so we can prove the following 
relation in (E^), 



i-i 



(ca + Ebjedi + {(-a + eb2)ed2 H h (fa + ef)„_i)ed„_i = '^{<^c, + Sd, ) ed^ 

3=2 k=l 



(43) 
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That means the matrices in the above relation is still in the graded tensor product form which agrees 
with the definition in (|2^). And the above defined nested transfer matrix can be defined as the supertrace 
on the auxiliary space for the reduced monodromy matrix which satisfies the Yang-Baxter relation. 

t«(u,K})=strT(i)(u,K}) 

R't^{u-v)T['\u,W})T!,'\v,W})=T^'\vAvu})T['\u,{v^})R^ (45) 

Thus we deal with almost the same problem as the original one. Repeating the graded algebraic Bethe 
ansatz method with grading FFB, we can find 

Ad) (.) = n ( ^^^-^ ) n - -'^^ + ^) + n ( ^ ) n - ) ^^^h-) • (46) 



j = l ^ ' k=l j=l ^ ' ■' ' k=l 

The eigenvalue A(^)(u) for grading FF can be obtained almost in the same way for the normal six- vertex 
model. We have 

n 0^ n (« - « -) - n ^ n (» - «>■ («> 

j = l ^ ^ ' k=l j = l ^ ■' ' k=l 

In addition to these relations, we denote by A(u) the eigenvalue of the original transfer matrix t(u). We 
have 

A(^) = n ( ^-^-^ ) + *)^ + n ( ) """^^'^ (^)' (48) 



where Vj^^j^Vj should satisfy their corresponding Bethe ansatz equations. To summarize, the Bethe 
ansatz equations are listed as follows, 

Vk^\ ^ TTM^^^^fe^ T-r Vk-Vl-^ 

m i m^^^ . m'^^ i 

Mfc - ^'j - 5 7=1 -^^J-^ n-^^k~i' 

n '^'^^'f = n (51) 

where k take values from 1 to m, m*^^) and m^^' in (]49|) , (pO|) , (|5l|) respectively. Note that we have already 
redefined the parameters Vj,^j,Vj by shifting and ^ respectively in the above relations. According 
to the definition of Hamiltonian in section 3, we can find the energy of the Hamiltonian H''^ as 

m -. 

E = N-Y^-^. (52) 

J=l J ^ 4 
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We shall use the notations introduced in ref. |Tj . We denote respectively iVf and Ni as the number of 
single electrons with spin up and spin down, and Ni, number of local electron pairs. As we have pointed 
out, the reference state is chosen to be | ti>: ^^^d the vacuum state is the full-filled state. Therefore, 
an electron means a hole. The numbers m, m^-^^ , m*-^) appeared in above relations can be written as 
m = N-^ + Ni + Ni, TO^^^ = iV| + Ni, m*^^^ = iVj^. Since the Hamiltonian commutes with the Hubbard 
term Xj, this term will not change the eigenvectors and the Bcthe ansatz equations. 

5 Analysis of the Bethe Ansatz Equations and Discussions 

Using the standard method [l2| , we next analyze briefly the Bethe ansatz equations. By taking the 
logarithm of the Bethe ansatz equations, we obtain the following set of equations, 

NHvk, 2) - I] - VI, l)-Y. "^^f^i - ^fe' 2^ = ^^^^ 
1=1 j=i 

J2 -VI,-)- ~ ^) ~ ^('''^ " 2^ " ^^^^ 

;=1 /3=1 a=\ 

^ $(^^ - vb, " 1) = 27ri^6, (55) 

/3=1 c=l 

where $(A, a) = 2arctan{X/a). Ik,Ja Siiid Ki, are integer or half-odd integer depending on m^nS^^ and 
m^'^\ In the thermodynamic limit N — > 00, Vj,^j,Vj become continuous varibles, Ik,Ja and Kb are 
connected with the distribution functions of roots and holes, the Bethe ansatz equations are rewritten in 
the following form 

<^{v) + (Jh{v) = Ki/2{v) - / dv'Ki{v - v')a{v') + / d^Ki/2{v - ^)u){^), 

J-B J-BW 

uj{n) + LUhifJ-) = / dvKi/2{n - v)a{v) - / dfi' Ki{n - fi')uj{n') + / dvKi,2{li - v)t{v), 

J-B J-BW J-_B(2) 

T{iy) + Th{iy) ^ dfiKi/2{iy ~ fiMfi) - dv' Ki{v ^ v')t{v'), (56) 

where Ka{\) = ^(x"+a'i) ■ These results are completely the same as those for SU(4) spin chain, see ref. pO| . 
The energy is written as 

E = N-2tt[ Ki/2{v)a{v)dv. (57) 

J-B 

The ground state energy per site is obtained as 

Eo = l-{hn2+^). (58) 

Because the coupled integral equations ( p6| ) are the same as for the case of normal SU(4) spin chain, the 
properties of the low-lying excitations of the system under consideration remain to be the same as those 
in ref.[f0|. 
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In conclusion, we have presented in this paper the Bethe ansatz equations for the SU(4) extension of 
the Hubbard model. We have proved that the Hamiltonian is equivalent to the graded SU(4) permutation 
operator. Using the graded algebraic Bethe ansatz method, we have found the eigenvalues of the transfer 
matrix. 

In this paper, we have only dealt with FFBB grading, and there are only one matrix representations 
for the Fermion and Boson generators. We can consider other representations as presented in ref . pl| and 
different gradings. The SU{4:) extension of the Hubbard model is different from the SU{2\2) case, but we 
use in this paper the same 16 generators as in ref. to represent the Hamiltonian. The structure of the 
Bethe ansatz eigenvectors for 5'C/(4) case seems to be the same as the SU{2\2) case. We can also give a 
detailed calculation of the thermodynamic Bethe ansatz for the SU{4) case as that of 5'£/(2|2)|]lT|. 
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